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The chameleon scalar field is a matter-coupled dark energy candidate whose nonlinear self- 
interaction partially screens its fifth force at laboratory scales. Nevertheless, small-scale experiments 
such as the torsion pendulum can provide powerful constraints on chameleon models. Here we de- 
velop a simple approximation for computing chameleon fifth forces in torsion pendulum experiments 
such as Eot-Wash. We show that our approximation agrees well with published constraints on the 
quartic chameleon, and we use it to extend these constraints to a much wider range of models. Fi- 
nally, we forecast the constraints which will result from the next-generation Eot-Wash experiment, 
and show that this experiment will exclude a wide range of quantum-stable models. 



I. INTRODUCTION 



Evidence for an accelerating cosmic expansion is now 
solid [3-13l- The simplest explanation for this accelera- 
tion, a "cosmological constant" vacuum energy density 
A7l/pj which does not interact with Standard Model par- 
ticles, remains consistent with the data; however, it must 
take an extremely small value ~ 10 — 120 Mpi which is dif- 
ficult to explain without fine-tuning. Alternative theo- 
ries (5l-[l5| explain the smallness of this number dynam- 
ically, either through tunneling to a low-energy vacuum 
or through a slow reduction of the vacuum energy known 
as "degravitation." Since the simplest of these models re- 
duce at low energies to effective scalar field theories pos- 
sibly coupled to known particles, it is interesting to con- 
sider such a scalar field "dark energy" independently of 
these more fundamental theories, and to ask how it may 
differ from a cosmological constant. Generally speak- 
ing, such differences take two forms: a slow evolution of 
the background energy density fill Il7j ; and couplings to 
Standard Model particles, which we consider here. 

Large fifth forces have not been observed, so coupled 
dark energy must have some mechanism to screen such 
couplings at laboratory and solar system scales. Galileon 
fields invoke the nonlinear Vainshtcin mechanism to re- 
duce their effective couplings at high densities [I3 - |20| . 
Symmctron models decouple from matter through a sym- 
metry restoration at high densities, while fifth forces exist 
in a symmetry-broken phase at low densities j2l| - [23j . The 
first screened scalars to be considered as dark energy can- 
didates are chameleon models, which evade constraints 
by becoming massive in high-density environments [24l - 
|26[. The current article focuses on chameleon models. 

Although these scalar fields are screened, such screen- 
ing mechanisms are not perfect. Laboratory experiments 
are powerful probes of residual fifth forces and new parti- 
cles which could result from coupled dark energy [23, Ull ■ 
Particles of a photon-coupled scalar could be produced 
through oscillation in a background magnetic field and 
detected using "afterglow" experiments (29l - [36j . Fifth 
forces may be probed directly through small-scale tests 
of gravity such as torsion pendulum experiments and 
Casimir force measurements I28L I37H451 . 



The goal of this article is to use torsion pendulum ex- 
periments such as Eot-Wash [3^| to constrain fifth forces 
from chameleon dark energy models [4(| I46l - l5l1 | . Pre- 
vious work [50j | used the numerical computations of \at\ 
to calculate the three-dimensional field configuration di- 
rectly for the geometry of the Eot-Wash experiment, a 
powerful probe of gravitation-strength fifth forces at sub- 
millimeter scales. In this work we develop a simple, ac- 
curate estimate of the field profile for such experiments 
by approximating the matter distribution locally as one- 
dimensional and planar. This one-dimensional plane- 
parallel (lDpp) approximation allows us to compute the 
field on the surface of the source and test masses in a tor- 
sion pendulum experiment, from which the energy and 
torque can be found. We show that the lDpp approxi- 
mation agrees with the numerical calculations of |47j and 
the data analysis of [5(| for Eot-Wash, and we estimate 
the constraints on a much wider range of models. 

Recently it was shown that a subset of chameleon 
models is "quantum-stable" in the sense of having small 
one-loop corrections to the effective mass and bulk field 



value |51| . For gravitation-strength couplings, quantum- 
stable models lie just beyond the bounds of Eot-Wash. 
Using our lDpp approximation, we forecast constraints 
from the next-generation Eot-Wash experiment and show 
that it rules out a large range of such models. 

The paper proceeds as follows. Section [TT] introduces 
the chameleon model and its fifth force screening mecha- 
nism. In Sec. IIIII we study in detail the one-dimensional 
planar problem, which is exactly solvable for power law 
chameleon potentials. Solutions of this one-dimensional 
problem are used to approximate the expected torsion 
pendulum signal in Sec. IIV1 and Sec. IVl concludes. 



II. CHAMELEON FIELDS 

A. Equation of motion 

The chameleon field <j> is a canonically normalized 
scalar field with a nonlinear self-interaction and a mat- 
ter coupling mH^. A simple matter interaction results 
from the conformal coupling of the chameleon field to the 
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metric, as given by the action 

{d<t>? 
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in the flat-spacctime case appropriate to laboratory tests. 
Here V((j>) is the chameleon potential, and matter is rep- 
resented as Fermion fields i^i with Lagrangian density 
C m . Conformal coupling of the chameleon results in a 
universal coupling constant /3 > to all Fcrmionic mat- 
ter, a feature which is stable with respect to quantum 
corrections [52ll53l|. In a background matter density p(x) 
the chameleon equation of motion is 



dVcs 
d<f> 



Mi 



-p(x)<p, 



pi 
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where V e g is the effective potential. We have neglected 
terms of higher order in /3</>/Mpi since this quantity will 
be small in all models of interest. 

In this work we will primarily be concerned with the 
static case, in which the equation of motion reduces to 



V 2 <P = V'(<j))+ Pp/Mi 



pi- 



(4) 



Deep inside a bulk medium of constant density po even 
the spatial derivatives vanish. The bulk field 4>b(po) then 
satisfies V ((/>b(po)) + Ppo/Mpi = 0. The effective mass 
associated with small fluctuations about a field <p is given 
by m cff (0) 2 = V" {<!>). 

The chameleon potential V{<f>) must be chosen to fit 
the cosmological data and to screen fifth forces locally. 
Cosmology requires that V > vary sufficiently slowly 
with time, and we will see that the chameleon effect re- 
quires V < and V" > 0. Since the cosmic acceleration 
is sourced by constant or slowly- varying parts of V while 
laboratory experiments are sensitive only to derivatives of 
V, we choose a constant-plus-power-law potential which 
splits these two regimes: 



(5) 



V{<t>) = Ml 1 + 7 



Here Ma = 2.4 x 10~ 3 cV is the dark energy scale, so 
that the constant term M^ drives the cosmic accelera- 
tion. The second term, in which 7 > and n are dimen- 
sionlcss numbers, can be probed in laboratory experi- 
ments. For n = 4 it is conventional to define A = 4!7. 
The bulk field and mass are given by 



4>b(p) = cr n M A 



rip 
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m eS (p) = M A |n-l|2(|n|7)a«-s 



where a n = sign(l — n). 
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This potential is the large-field limit of the exponential 
potential V = M% exp(i<f> n /M]t) of [2(j. Henceforth we 
work with ([5]) whenever specific examples or constraints 
are presented. Furthermore, the mass does not increase 
with density when 1 < n < 2, and we will see that ([5]) 
is constrained by cosmology when —1/2 < n < 1. Thus 
we only consider models with n < —1/2 or n > 2. Note 
that, due to our sign convention j3 > 0, the field <j) w iH 
be negative for n > 2 and positive for n < 0. For all such 
n, 4>b(p) decreases as p increases. 



B. Chameleon and thin-shell effects 

Chameleon phenomenology is characterized by the 
presence of two regimes: a linear, or "unscreened," 
regime; and a nonlinear, "screened" regime. In the linear 
regime, the potential derivative term on the right hand 
side of (@| is negligible, so the equation of motion is ap- 
proximately linear in </>. Furthermore, the source term 
remaining on the right hand side is proportional to p; 
thus (|1]) in the linear regime is similar to the Poisson 
equation V 2 * = p/(2Mp 1 ) for the gravitational poten- 
tial 'J. Since gradients of </> and $ vanish far from an 
object, 4> is equal to 2/3Mpi\f' up to an additive constant, 



A0 (lin) = 2/3Mpi A<J 



(8) 



where the A denotes a difference between two spatial 
positions. The linear regime applies, for example, to a 
dense object of sufficiently small volume in a sufficiently 
low-density bulk. 

Suppose that the volume of such an object is increased 
at constant density. For a characteristic size r the grav- 
itational potential ^> oc pr 2 , and cf> will change with ^> 
throughout the linear regime. As 4> changes from its min- 
imum, V'((j>) will become large and negative, partially 
cancelling the source density on the right side of 

This cancellation, known as "screening" of the source, 
is characteristic of the nonlinear regime of chameleon 
models. In the nonlinear limit this screening becomes 
complete and the field asymptotically approaches its bulk 
value 4>b{p)- Since the gravitational potential continues 
to grow linearly, in the nonlinear regime the change in 
the field value will be much smaller than the linear ap- 
proximation ((SJ, 



Acj) 



(ni) 



< 2/3M P i I A* I 



(9) 



Suppose that the chameleon field at the center of an 
object in the nonlinear regime is </>o, and the field far away 
is 0oo. The gravitational potential at the center is ^q; far 
away, W is defined to be zero. For chameleon models with 
negative n, 4> ps 0, so © becomes \4>oo\ <C 2/?Mpi|\I>o|- 
For n > 2, the opposite is true; 0o ^> 0oo, so (O implies 
1 0o I •C 2/9Mpi|\Po|- The gravitational potential of a disk 
of radius rdisk and thickness Zdisk is approximately ^ ~ 
P?" d i sk z d i sk ; Ref. 53 finds * = pr disk z disk / ftM^) to be a 
good approximation. 
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FIG. 1: Mass scalings of various chameleon models. Labo- 
ratory constraints are most suitable for probing models with 
n < —1/2 and n > 2, while cosmological fifth force con- 
straints probe the models — l/2<n<lin which m e ff grows 
most rapidly with density. 



The chameleon fifth force in the nonlinear regime is 
suppressed by two effects known as the chameleon and 
thin-shell effects. The chameleon effect is the rapid 
growth of the effective mass m e ff with the size and density 
of a source object, which effectively converts a long-range 
force into a short-range one. For example, the Compton 
wavelength of a 7 = 1, n = — 1, /3 = 1 model increases 
from ~ 100 pc at cosmological densities to ~ 0.1 mm 
at laboratory densities. The thin-shell effect is due to 
the screening of the interior of a source mass. If the 
chameleon field is near its bulk value inside an object, 
and changes only in a thin shell of matter at the outer 
edge of that object, then a test mass outside that object 
will "see" only the fifth force due to that thin shell of 
matter. Due to the thin-shell effect, a chameleon with 
a gravitation-strength coupling (3 ~ 1 can easily be con- 
sistent with solar system fifth force constraints. Because 
of the chameleon and thin-shell effects, a model with ef- 
fective mass TO e ff at a given density is best probed using 
objects of size ~ m~g separated by a distance ~ m eg ■ 

Using the mass scaling (0, m e ft oc p(™~ 2 )/( 2 ™~ 2 ) 5 
we can estimate whether a given chameleon model can 
best be constrained by laboratory or cosmological ex- 
periments. Small-scale tests of gravity can exclude un- 
screened gravitation-strength fifth forces on length scales 
> 1 mm at densities p ~ 1 g/cm 3 , while cosmolog- 
ical probes exclude such forces on megaparsec length 
scales at the cosmic background density ~ 10~ 30 g/cm 3 . 
These approximate excluded regions are shaded in Fig- 
ure [TJ Suppose we have a model which can barely be 
probed in the laboratory, m c g ~ (1 mm) -1 , such as the 
models with n = 4, n = — 1, and n = —4 shown in 
the figure. This model will be too massive to probe 
cosmologically if the mass at cosmological densities is 
greater than (1 Mpc) -1 . This condition is satisfied if 

7^<TS lo Sio Tmin' implying either n < -1/2 or n > 1. 
The remaining models have rapid mass scalings and are 
best probed cosmologically. 



C. Quantum stability condition 

Reference (5lj derived a condition for the quantum sta- 
bility of a chameleon potential, that is, the condition 
that quantum corrections to the potential be subdom- 
inant over the range of field values probed by a fifth 
force experiment. Large masses, which help a chameleon 
model to evade fifth force constraints, also lead to large 
quantum corrections. For a general class of potentials, 
Ref. [5l| found the quantum stability condition m e g(p) < 
0.0073(/3/10 g/cm 3 ) 1 / 3 eV. Quantum-stable models are 
an interesting subset of all chameleon models, and we 
will see that the next-generation Eot-Wash experiment 
can exclude a range of quantum-stable chameleons. 

For a potential such as (O , the quantum stability con- 
dition is the requirement that the one-loop Coleman- 
Weinberg corrections to the slope of the potential and to 
the chameleon effective mass (that is to V and V") be 
smaller in magnitude than their tree-level counterparts. 
The one-loop correction to the potential is 
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log 
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where primes denote derivatives of the potential ([5]) with 
respect to 0, and po is a mass scale which characterizes 
the chameleon field in the experiment. Corrections to 
V and V" are the first and second derivatives of (QI 
respectively, so the quantum stability conditions are 



AV{_ 



loop 



V 



AV('_ 



loop 



V" 



< 1. 



Ill) 



Although we can choose po to make Ayi_i oop zero at 
any given field value, a fifth force experiment will probe 
a range of field values. <f> and m e ff can vary by an order 
of magnitude or more in a typical experiment, so quan- 
tum corrections will not be zero everywhere. When we 
consider the quantum stability of a model in a specific 
experimental setup in Section IPVl we will choose po from 
among the chameleon masses in the experiment so as to 
minimize quantum corrections. 

As an estimate of quantum corrections, we may set the 
log term in (fTOj) to unity and evaluate (fTTj) at the bulk 
field 4>b(po) corresponding to some density po- For the 
potential ([5]) the stability conditions are 
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\n(n- l) 2 (n-2)| 

327T 2 

|n(n- 2)(2n- 5)| 



(12) 



(13) 



Note that quantum stability imposes an upper bound 
on the self-coupling 7 for n > 2 and a lower bound for 
n < 0. In the case n = 4 the density-dependent term 
disappears, and (|12j|13p imply 4! • 7 = A < 32tt 2 /3 ps 
105. Casimir force constraints can rule out quantum- 
stable n = 4 chameleons with strong matter couplings 
/3 > 10 6 [ioj ]; however, quantum-stable models with 
gravitation-strength couplings /3 ~ 1 remain allowed. 
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III. FIELD PROFILE IN PLANAR SYSTEMS 



B. Planar gap 



A. Planar slab in vacuum 

The field profile in the vacuum outside an infinitely 
thick planar slab can be found exactly [4(| S3- Let 
p(z) — po<d(—z) where is the step function. Thus p 
is positive for z < and zero for z > 0; the face of the 
slab is the xy plane, and its normal is z. The static equa- 
tion of motion (j4|) in the vacuum, cPcft/dz 2 = dV/d<f>, is 
solved for the potential §5§ by 



z) = SV 1 + d -(n - 2) 2 7 M|-™0 SV 2 z 



(14) 



where 



is the field value on the surface z = of the 



slab in vacuum, 
equation of motion, we find 



Using 



2 dcf> dz 



to integrate the 



1 fd<j) 

2 Ydz 



Vc 



efiWsv 



,p) - V ef£ (4>i,p). 



(15) 



Choosing 0.; = 0b (po), corresponding to z — > — oo and 
P = Po, yields one equation relating SV to d(p(0)/dz; 
choosing 0; = 0b (0), corresponding to z — s- oo and p = 
0, yields another. Combining the two, and noting that 
d<p/dz —t as z —> ±oo, gives the simple result 



1 ) 0b (po)- 

n ' 



(16) 



An exact, closed- form solution is not available inside 
the thick slab. However, we can linearize the equation of 
motion around 0b (po) and require <fi(z) to be continuous 
at z = 0: 

0thick(z) « 0b(po) + [0 S v - 0B(po)]e mrff( " o)z . (17) 

The case of a thin slab is slightly more complicated. 
Suppose that the slab is centered at z = z c with a 
half-thickness of 5z. Guess a value 0c = (j>(z c ). The 
equation of motion linearized about 0c , and its solution 
<50thin(z) = 0(z) — 0c inside the slab, arc 



dz 2 

<50thin(z) 



V c ft,<f,(<l)c, Po) + m cS ((f>c) 2 
V e n,<t>((pG, Po) 



m e ff(0c) 2 



cosh(m e ff(0c)-2) - 1] 



(18) 
(19) 



Replacing SV by 0c + S(j) t h in (z c + Sz) in (fI3]l. we can 
find the exterior solution corresponding to this guess 0c • 
When the correct value of 0c is chosen, d(j>/dz as well as 
will be continuous at z c + Sz. However, if we consider 
thicker and thicker slabs, we cannot be sure that such a 
solution will exist. 



Consider a planar gap with p = p v bounded on the left, 
z < 0, by a thick slab of density p m L, and on the right, 
z > Az, by a thick slab of density p mR . That is, p(z) = 
PmL@{-z) + p v @(z)@(Az - z) + p mR 9(z - Az). There 
are four unknowns: the surface field values 0l and 0r at 
2 = and Az, respectively; the maximum field value g 
inside the gap; and z g , the point at which = g . 

Equation [15] can be applied to any interval 
[0(zi), 4>(zj)] over which p is constant. The intervals 
[<h(Pmh), 0l] and [0L,0 g ] give, respectively, 

2 



dz 



V(fo)-V(MPmL)) 
. /3PmL(0L - 0B(PmL)) 



M, 



pi 



F(0 g ) - V(<p L ) 



/3p v (0 g - L ) 



AI- 



■*pi 



(20) 
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Adding the two yields a relation between 0l and g . 
A similar procedure can be applied to the plane on the 
right. Thus we have 



0L 



PmL0B(PmL) - Pv0 g . ^(0B (Pmh)) ~ V((p s ) 



PmL — Pv 

PmR0B (PmR) - Pv0 g 

R — 

PmR — 

Next, we apply (fT5|) to 



/3M P1 1 (/j mL - Pv) 

V(0 B (PmR)) " U(0 g ) 



/3M pi (p mR - p v ) 



(22) 
.(23) 



for some arbitrary z 
between and z g in order to find d<p/dz inside the gap: 

g = V2 ] jv^)-V(^) + ^-^~cj )g ). (24) 

We can integrate with respect to between 0l and g 
to determine z g in terms of g : 



d^j^Mt 
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(-iV^r(i)inl 1 



=U 3 = 
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2 gap (0 g ,0L,Pv) 



<Mp v ) 



i(ra-l) 



xg. 



l '2 +l 



(25) 



where B x (a,b) = f Q t a x (l — t) b is the incomplete 
Beta function, and we have defined the shorthand Z aar> 



5 



for this expression as a function of the gap field <fi g , sur- 
face field </>l, and gap density p v . Repeating this proce- 
dure for the right side of the gap, z g < z < Az, 

Az- z g = Z gap (</> g ,0 R ,p v ), (26) 

which is similar to (|25p but with </>l replaced by </>r in 
the incomplete Beta function. 

Thus we have 0l , 4>r , and z g in terms of <p g , while the 
sum of ([2"5"]) and ([2"B]) implicitly defines <fi g as a function of 
the known gap size Az. By guessing </> g , comparing the 
resulting Az to the known value, and itcratively refining 
our guess, we can find (p g . Once </> g is known, we can 
integrate (|24|) to find </>(z) within the gap. For example, 
given any 0o = <fi( z o) between 0l and </> g , we obtain for 
z g — zq a formula similar to ([25]) with 0l replaced by 4>q. 

The series sums in (|25l [26]) will converge quickly un- 
less 4> g is close to </>b(/9 v )- In that case, the fifth force 
on each slab will be small anyway; the gap is large 
enough that the opposite slab does not pull </> g very far 
from its bulk vacuum value. When fifth forces are large, 
even the i = term alone is a reasonable approxima- 

tion: ^m^WAz « B 1 _^ M n(l/2, 1/2 - 1/n) + 

£>!_(0 g /0 R )n (1/2, 1/2— 1/n). Henceforth, in our numeri- 
cal calculations, we truncate 2 gap after the i = 5 terms. 

To summarize, we have shown how to compute the sur- 
face field </>r(Az) as a function of gap size Az in a planar 
system. This result will be essential to our approximation 
for torsion pendulum experiments in Sec. IIVI 



C. Thin slab in planar gap 

Before proceeding to experiments we study one final 
planar configuration, the thin slab inside a planar gap. 
This will allow us to calculate the chameleon screening 
caused by the electrostatic shielding foil between source 
and test masses in short-range fifth force experiments. 

Figure [5] shows the geometry considered here. A cen- 
tral slab of width Azc and density p m c sits between two 
thick slabs, one on the left at a distance Azl with a den- 
sity /9 m L, and another on the right at a distance Azr 
with a density p m R. Gap 1, between the left and cen- 
tral slabs, has a "vacuum" with density p v i, while gap 
2, between the central and right slabs, has a densty p V 2] 
we assume p v i,Pv2 -C p m L,PmC,PmR- Field values deep 
inside the left and right slabs are 0b l = (psiPmh) and 
4>bk = 0B(PmR,), respectively. Without loss of generality, 
let z = be the face of the left slab. Given these values, 
we wish to find: 0li = 0(0); g i, the maximum inside 
gap 1; <f>ci = 0(Azl); 0c 5 the minimum inside the cen- 
tral slab; <pc2 = 0(Azl + Azc); 0g2, the maximum inside 
gap 2; </>R2 = 0(Azl + Azc + Azr); and the positions 
z g i , z c , and z g 2 at which the local extrema (j) g i , </>c , and 
</> g 2, respectively, are attained. 

For these ten unknowns we have ten equations. Two 
are obtained by applying (|15l) to itcrvals [0Li;0gi] and 




FIG. 2: Density (shaded yellow region) and field profile (solid 
blue line) for a thin planar slab inside a gap between two thick 
planar slabs. The horizontal axis shows the z coordinate, the 
distance normal to the planes, while the vertical axis shows p 
and <j> in arbitrary units. Distances, densities, and field values 
are labeled. 



[4>gi, 4>Ci\] two more by evaluating the thin-slab lineariza- 
tion (|19p at z = Azl and Azl + Azc; and two more by 
applying ([15]) to itervals [0c2,0 g 2] and [0 g 2,0R2]- The 
final four are found by applying (|251 [26]) to the intervals 
[0, z g i], [z g i, Az L ], [Az L + Az c , Az L + Az c + z g2 ], and 
[Azl + Azc + z g 2, Azl + Azc + Azr], that is, to the left 
and right sides of gaps 1 and 2. 



= 


= Kff(< 


^PrnL)\%l + V eS (<t>,P V l)\Z 


(27) 


= 


= Kff(< 


^P V l)\Z + V eS (<t>,PmC)\i C cl 


(28) 




4>G1 = 


(j> C + <5</>thm(Az L ) 


(29) 




4>C2 = 


0c + <50thin(Az L + Az c ) 


(30) 


= 


■■ v cS {< 


^ Pm c)\tc+VMP V 2)\% 2 2 


(31) 


= 






(32) 


Z g l = 




^gl,0Ll,Pvl) 


(33) 


Az L = 


Zgap(<, 


^>gl,0Cl,Pvl) + %1 


(34) 


Z g2 = 


Zgap(<, 


4g2, 0C2, Pv2) + Az L + Az c 


(35) 


Az R = 




^g2, 0R2, Pv2) + Z g2 - Az L - Az C 


(36) 



Note that if <fi g i, g 2, and 0c are specified, then ([2~?1 
I36|) immediately give the other seven unknowns. We also 
obtain gap sizes Az^ and Az R and slab thickness Az c ; 
however, these will not necessarily match the givens Azl, 
Azr, and Azc- In order to find the correct g i, g 2, and 
4>c , we minimize (Az^ — Azl) 2 + (Az R — Azr) 2 + (Az' c — 
Azc) 2 numerically. 

In the symmetric case, p m L = p m R ; Pvi = Pv2, 
and Azl = Azr, the problem simplifies consider- 
ably. Matching the field derivative at the surface of 
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+ thin-slab linearization 
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• numerical integration (CLN) 
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m eff Az foil 

FIG. 3: Factor / aU p by which the chameleon fifth force is sup- 
pressed by a shielding foil. For a large range of chameleon pa- 
rameters 7, n, and /3, the suppression factor, shown by points 
on the plot, is a function of the foil thickness A.Z{ n in Comp- 
ton wavelengths. / sup is well-approximated by the function 
sech(2m c ff Zf ii) (solid line). The thin-slab calculation of / sup , 
from (|27I36[) , is shown by blue "+" -shaped points; the scatter 
at large m e ffA2f ii is due to numerical error, since A(/i(foil) is 
the difference of two nearly equal numbers. Solid black circles 
show / SU p computed directly from arbitrary-precision numer- 
ical integration of the equations of motion. 



at its center will be close to its bulk value, and the field 
on the side facing the test mass will be very weakly de- 
pendent on the field at the opposite face. In the case of a 
thin slab, however, this screening will not be total. The 
results of Sec. IIII Cl provide an excellent approximation to 
the surface field in the presence of a shielding foil. The 
change in source mass position will result in a change 
A0(foil) in the field on the surface of the test mass. 

Figure [3] shows the suppression factor / sup = 
A0(foil)/A0(no-foil) for a Azq = 10 fim foil at the center 
of a gap with Azl + Azr + Azc = 100 /im. The density 
of each slab is 10 g/cm 3 and the density in the gaps is 
10 _12 g/cm 3 , corresponding to air at room temperature 
and a pressure of < 10 -6 torr. A</>(foil) is found by vary- 
ing Azl by 10 ^m in either direction and using (j27H36[) 
to find the change in 4>R2 ■ /sup is approximately equal to 
sech(2m c ff Azc), where m e ff is evaluated at </>B(Pmc)- 

At large m e ffAzf ii, / sup is the difference between two 
nearly-equal numbers. The resulting numerical error is 
responsible for the scatter in the "+" -shaped points at 
TOeffAzfoii > 5 in Fig. [3J In order to verify our ap- 
proximation for / sup in this regime, we integrated the 
equation of motion numerically using the CLN arbitrary- 
precision arithmetic package [54[ . The resulting / sup val- 
ues, shown as filled circles in the figure, agree well with 
/sup ~ sech(2m e ffAzf oi i). 



the central slab gives (f5^ hin (Azc/2)) 2 /2 = V e s{4>C + 
<50thin(Az c /2),p v i) - V e fi(<j>g,i,Pvi), where the prime (') 
denotes d/dz. Thus (f> g i determines <f>c. We need only 
solve numerically for the value of <fi g i which gives the 
right gap size Azl using (|3"3l IM)) . 

D. Force suppression due to shielding foil 

Forces at short ranges between small source and test 
masses in a fifth force experiment will typically be dom- 
inated by electrostatic effects. Even electrostatic forces 
between stray charges can swamp gravitation-strength 
forces. Thus most such experiments stretch a thin, 
grounded conducting foil between the source and test 
masses to shield the test mass from these electrostatic 
forces. Such a foil can screen the chameleon fifth force as 
well, so we study it here. 

First, consider a system with two thick planar slabs. 
Let the slab on the left be the source mass and the one 
on the right the test mass. A change in the position of 
the source will change the surface value of the field on the 
test mass by some amount A0(no-foil). The change in 
the fifth force on the test mass, the experimental signal, 
is proportional to A0(no-foil). 

Now suppose that another slab, corresponding to the 
shielding foil, is placed between the source and test 
masses. In the limit that this central slab is thick, it 
will completely screen chameleon fifth forces. The field 



IV. TORSION PENDULUM EXPERIMENTS 
A. 1-D plane-parallel approximation 

Thus far we have studied planar configurations because 
exact solutions exist. However, torsion pendulum exper- 
iments such as Eot-Wash do not measure fifth forces in 
the i direction between planar slabs. A better approx- 
imation is a pair of slabs with features such as holes or 
grooves. The z positions of both slabs are fixed, and the 
source slab is moved in a direction parallel to the planes 
which we call x. As features in the source mass move 
past those in the test mass, forces are exerted in the ±x 
direction. 

Reference [13] computed the field numerically for such 
a density configuration. The field <p was discretized on a 
three-dimensional grid of points and the Hamiltonian was 
minimized with respect to this discrete set of field values. 
Figure 3] shows the field configuration when the features 
on the source and test masses are circular holes. Once 
the field is known, the force on the test mass, occupying 
a volume V, can be computed directly from the gradient 
of the interaction potential, F = — J v d 3 x(V(f))pMp ] L . 
However, solving for the field over a range of x positions, 
for large ranges of 7, n, and /3 values, is computationally 
expensive. Accurately accounting for force suppression 
due to the shielding foil requires discretizing space on 
length scales <C 10 urn, yet covering a region of size ~ 
10 mm. Thus Ref. [50[ restricted itself to n = 4 and 
/?< I- 
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FIG. 4: Chameleon field for a hole in the source mass moving 
past another one in the test mass. The chameleon potential 
is V(tf)) = jf</> 4 with A = 1, equivalent to ([5| with n = 4 and 
7 = 1/24; the matter coupling is /3 — 1. 



Here we make a series of approximations to which we 
refer collectively as the one-dimensional plane-parallel 
(lDpp) approximation. Consider two parallel slabs as 
in Figure [5] (Top), with a hole or groove in the lower 
(source) slab. Our goal is to estimate the field at each 
point on the surface of the source and test masses by ap- 
proximating the matter distribution near that point as a 
planar gap. We use the results of Sec. IIII Bl which found 
the surface field </>r (Az) as a function of planar gap size 
Az. At any point on the surface of either slab, let Az 
be the distance to the nearest point on the opposite slab, 
and approxiate the surface field as 0r(Az). Using this 
field profile, we compute the energy. The force in the x 
direction is the rate of change of this energy as a hole on 
the source mass passes by one on the test mass. 

Let the distance between the source and test slabs be 
Azg_T, and consider points p, q, and r on the surface 
of the test slab, as shown in Fig. [S] (Top). Our one- 
dimensional plane-parallel calculation makes the follow- 
ing approximations. 

1. Each slab is thick enough that the chameleon at- 
tains its bulk value deep inside. 

2. The field at points such as p, which are not directly 
across from the hole on the opposite slab, is equal to 
the surface field </>r in a one-dimensional planar gap 
of width Az(p) = Azs-t as studied in Sec. IIII Al 

3. The field at a point such as q or r, which is directly 
across from the hole in the opposite mass, is equal 
to the surface field 4>R m a planar gap with Az 
equal to the distance to the nearest point on the 
opposite slab. For example, <j> at q is equal to the 
surface field in a gap of size Az(q) shown in the 
figure. 

4. Since the surface field should not change on length 
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FIG. 5: One-dimensional plane-parallel approximation. 
(Top) A feature, such as a groove or hole, on the source mass 
faces the test mass. Three points, p, q, and r, are labeled on 
the test mass, and Az(p), Az(q), and Az(r) are respectively, 
the distances between each of p, q, and r and the nearest point 
on the source slab. (Bottom) lDpp approximation at r. 



scales larger than the Compton wavelength, we ne- 
glect the energy due to a transition region of width 
m~g at the edge of the region across from a hole. 

5. When computing the total energy of a configura- 
tion, only the field inside the source and test masses 
will be counted; changes in the field profile inside 
holes in each disk, as well as in the gap between 
disks, arc neglected. 

6. A shielding foil of thickness Azi n between source 
and test masses reduces the force by a factor of 
/sup = sech(2m c ff Azfoii) where m e ff is the effective 
mass of the chameleon field at the bulk density of 
the foil. 

Figure [5] (Bottom) shows this lDpp approximation at 
point r. The geometry of Fig. [5] (Top) at r is replaced 
by a one-dimensional planar gap in which the field can 
be calculated simply. 

Several of these approximations cause us to underesti- 
mate the force somewhat. In particular, the approxima- 
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tion [3] above adds matter near the opposite mass. This 
means that even in non-overlapping regions, the field in 
the lDpp approximation will be closer to its bulk value. 
Thus the energy difference as the source mass moves is 
underestimated, leading to an underestimated force. Fur- 
thermore, aligning holes in the source and test masses will 
lower the energy associated with the field inside the holes 
and gaps as well as inside the material of the slabs them- 
selves. Approximation[5]ignores this energy change, lead- 
ing to an underestimate of the force. Approximations [T] 
and [5] also lead to slight underestimates. 

By its nature, the one-dimensional plane-parallel ap- 
proximation will predict no torque in Eot-Wash due to a 
massless field, hence no sensitivity to a 1/r 2 force such 
as Newtonian gravity. At large Azs_t we can use (jT4j) 
to approximate m e s ~ 0("~ 2 )/ 2 ~ 1/Azs-t- Thus the 
chameleon becomes effectively massless in the limit that 
Azg_x is much larger than the sizes of the features in the 
disks. Since Eot-Wash is sensitive to Newtonian torques, 
the lDpp approximation will underestimate the torque 
in this limit. Such an underestimate is not significant in 
Eot-Wash, whose chameleon constraints are dominated 
by separation distances Azg_x much smaller than the di- 
ameters of the holes. However, it does mean that the 
lDpp approximation will substantially underestimate the 
signal in a "chameleon lightning-rod" experiment such as 
that suggested by [55| . 

On the other hand, approximation [2] potentially leads 
to an underestimate of the energy at points such as p, 
which are not directly across from the hole on the source 
slab. If p is within a few Compton wavelengths of the 
edge of the hole, then the field there will be somewhat 
larger than expected for a gap of width Azg_x, hence 
its energy will be somewhat greater. Thus approxima- 
tion[5]incrcascs the energy difference between overlapping 
and non-overlapping regions, leading to an overestimate 
of the force. We shall see in Sections IIV CIIIV Dl that 
this overestimate is small for the current-generation Eot- 
Wash experiment but nontrivial for the next-generation 
experiment. Additionally, the size of the transition region 
in approximation U is just an estimate; it could be 1.5 
or 2 Compton wavelengths rather than one. Our choice 
above will lead to a slight overestimate of the force for 
the lowest chameleon couplings. 



B. Torsion pendulum 

Here we apply the lDpp approximation of the previous 
section to a hypothetical torsion pendulum similar to the 
Eot-Wash experiment (39j |. Such an experiment consists 
of a pair of parallel, rotating disks with matching holes 
at regular intervals, as in Figure [51 The lower disk, the 
source mass or "attractor," is mounted on a turntable 
which keeps it rotating uniformly. The upper disk, the 
test mass or "detector," is a torsion pendulum allowed to 
rotate freely. If there is a fifth force, then the test mass 
will experience small torques as the holes on the source 




FIG. 6: Geometry of disks used in Eot-Wash torsion pen- 
dulum experiment (from [47|; not to scale). (Left) Cur- 
rent experiment, with two rows of 21 holes in each disk. 
(Right) Next-generation experiment, with 120 radial grooves 
in each disk. 



mass move in and out of alignment with those on the test 
mass. 

The chameleon fifth force which results from moving a 
hole on the source disk past one on the test disk is the 
energy cost per unit distance of the change in the field 
configuration. Using the lDpp approximation, we can 
estimate the field configuration on the surface of each 
disk as one hole of radius rsh on the source disk rotates 
past another of radius rxh on the test disk. Assume that 
both disks have the same density p m . Let the origin of 
the coordinate system be the point on the source disk 
directly across from the center of the test disk hole, with 
z parallel to the rotation axis of the disks and x in the 
direction of motion of the source hole (that is, the tan- 
gential direction) . If the source hole is far away, then the 
field at a position r = \J x 2 + y 2 on the surface z = 
of the s ource disk is the surfa ce field in a gap of size 
Az = ^(Az s . T ) 2 + (r T h - r) 2 . Since Az > Az s _ T , the 
field </> will be greater on the portion of the source disk 
across from the test disk hole. Thus will be farther 
away from its energy- minimizing value </>b(Piii)- On the 
other hand, if the source and test holes overlap, then 
Az = Azs-T over the maximum possible area on both 
disks, minimizing the energy. 

Figure [7] shows the field profile <f>{x, 0, 0) on the surface 
of the source disk when the source hole is displaced from 
the test hole. The lDpp approximation is in close agree- 
ment with the three-dimensional numerical calculation 
of [l?) except for a transition region at the edge of the 
hole. Figure [8] compares the lDpp approximation and 
the 3D numerical calculation for a range of geometries 
and models. In all cases the two agree reasonably well. 

Now that the lDpp approximation has given us the 
field c/> B (x, y, 0) on the surface of the source disk, we may 
find the energy. Assuming that the disk is a thick slab, 
we approximate the field inside it using the thick-slab 
linearization (|17[) . Then the energy inside the region of 
the source disk across from the test mass hole, assuming 
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FIG. 7: Field on the surface of the source mass, with source 
and test holes offset, for V{4>) = jr</> 4 , A = 1, /J = 1, and 
Azs-t = 0.2 mm. The horizontal axis shows the tangential 
direction, and x = coincides with the center of the test mass 
hole. The lDpp approximation agrees quite well with the 3D 
numerical calculation of 14711. 





3 -4 -2 2 4 
x [mm] at surface al source disk 




(d) 


1 Dpp approximation 




-6 -4 -2 


2 4 
of source disk 



FIG. 8: Similar to Fig. [7] but with (a) a greater source hole 
displacement, (b) Azs-t = 0.1 mm, (c) A = j3 = 0.1, and 
(d)n = -l, 7=1, 13 = 1. 



that the source hole is far away, is 
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where </>b and m e g are evaluated at the disk density p m . 
The subscript h denotes the region of the source disk 
across from the hole. 

Next, consider a region of the source disk far from any 
hole on the test disk. In that region, the two disks will 
look like a pair of parallel planes with a separation Azs_t, 
so cf> s will be a constant on the surface. Let E n h be the 
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FIG. 9: Comparison between the lDpp approximation and 
the 3D numerical computation of [47| for a (f> 4 chameleon with 
A = ft = 1 and disk separation Azs-t = 0.1 mm. 



energy of a region of the same size, where the subscript 
nh is short for "no hole." Then AEc, = E^ — E n h is the 
energy cost in the source disk associated with each hole 
in the test disk. Similarly, we may compute A£t, the 
energy cost in the test disk. 

Finally, we may compute the total energy and torque. 
The amplitude of the total energy Etot will be half of the 
total energy change AE$ + AE^ , multiplied by the total 
number A^oies of holes, which is 42 for Eot-Wash. Let 9 
be the rotation angle, and define 9 = to be the angle 
at which source and test disk holes are perfectly overlap- 
ping. For equally-spaced holes in N TOWS = 2 rows, the 
frequency with which E to t varies is i>h = Nh i es /N I0WS . 
Multiplying by the force suppresion factor / SU p, we ob- 
tain the total energy and torque, 

-Etot = --N holes f sup (AE s + AE T ) cos(M) (38) 

t = ^hoies^/supCA^s + AEt) sm(v h e). (39) 

Figure [9] compares (l39l to the three-dimensional nu- 
merical computation of [47| . The lDpp approximation 
underestimates the torque by a factor of about 2.5. Fig- 
ure [10] demonstrates that this underestimate becomes 
worse by a factor of about two at separations around 
a few millimeters, the diameters of the source and test 
holes, as expected from Sec. IIV A[ However, since the 
torque itself falls off rapidly with separation distance, 
constraints will be dominated by small Azg_x- Thus this 
worsening of the lDpp approximation at large Azg_x will 
not have a significant effect on the final constraints. 



C. Constraints 

Using the lDpp approximation developed above, we 
may quickly estimate constraints on power law chameleon 
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FIG. 10: Comparison between the lDpp approximation and 
the 3D numerical computation of [13] for a cj> A chameleon with 
A = 1/10 and j3 = 1, as a function of Azs-t- 



models from the current-generation Eot-Wash experi- 
ment. This experiment has A^hoies = 42, N mws = 2, 
r-Th = 2.4 mm, and rgh = 1-6 mm. The source and 
test disks were made of molybdenum, with a density 
p m = 10 g/cm 3 , while the laboratory vacuum density 
was p v = 10 -6 torr ~ 10 -12 g/cm 3 . Eot-Wash probes 
torques over a range of disk separations; however, here we 
approximate the experiment as excluding torques greater 
than 0.003 fN-m at Az s _ T = 0.1 mm. 

Figure im shows our approximate lDpp Eot-Wash con- 
straints. In particular, Fig. [TTJ (Top) compares lDpp con- 
straints to the more precise numerical calculation of (5fj| 
for the 4 chameleon. In the range 0.01 < f3 < 1 cov- 
ered by both sets of constraints, the lDpp exclusion lower 
bound on A agrees well with the more precise calculation. 

The lDpp calculation underestimates constraints near 
P = 1 due to Approximation [5] of Sec. IIV A| which ne- 
glects the contribution to the total energy of the field 
in the gap between disks. We could potentially correct 
for this underestimate by including an extra factor in 
([551 r55|) and using the numerical computations of [13] to 
calibrate this factor. Meanwhile, around /3 = 0.01, the 
lDpp calculation overestmates constraints. Approxima- 
tion 0] of Sec. IIV Al assumes a transition region of width 
m~g associated with each test and source mass hole, but 
this is just an estimate. We could potentially include 
another factor parameterizing the number of Compton 
wavelengths in the transition region, and then adjust it 
to match Ref. |47| more closely. 

Since the goal of the present work is an estimate of Eot- 
Wash constraints rather than a rigorous data analysis, we 
do not fit these two "fudge factors" to [I?) ■ We have cho- 
sen an Eot-Wash exclusion limit of 0.003 fN-m such that 
the lDpp constraints approximately match those of [50| , 
which is equivalent to estimating a value for the first of 
these factors. Our choice is roughly consistent with the 
Eot-Wash data presented in (3t| , and the resulting lDpp 
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FIG. 11: lDpp excluded regions (shaded light green). The 
black shaded regions identify models which are linear in- 
side the source and test masses; these are excluded unless 
/3 is very small. (Top) Quartic chameleon, V(cf>) = jj<t> 4 - 
The long-dashed blue line shows the constraints of [3SL"l50l| . 
Models above the short-dashed purple line have large quan- 
tum corrections. (Bottom) Inverse power law chameleon, 
V{(j)) = M"a(1 + jMa/4 1 )- AH models shown have small quan- 
tum corrections. 



constraints are a slight underestimate in the strongly 
nonlinear regime /3 > 1. Meanwhile, we do not adjust 
the second factor at all. 

Models in Fig [TT] (Top) above the dashed purple line 
have large quantum corrections; they fail the quantum 
stability conditions discussed in Sec. Ill CI For a range of 
matter couplings 1 < /3 < 100, Eot-Wash excludes all 
quantum-stable n — 4 chameleon models. 

Constraints on the n = — 1 chameleon are shown in 
Fi gure [TT1 (Bottom). A more rigorous analysis such 
as [5 Of does not exist for this model. The quantum sta- 
bility conditions do not exclude any models shown here; 
quantum corrections are « 20% of the tree level values 
in the bottom right corner of the plot, and smaller else- 
where. For 7 = 1, the model in which the power law term 
in the potential ([5]) has the same energy scale as the dark 
energy, Eot-Wash excludes 0.01 < j3 < 15 in the lDpp 
approximation. 

Figure [15] shows constraints in the /3, n plane for several 
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FIG. 12: Excluded region in the /3, n plane, for 7=1 (green 
shaded region), 7 = 0.1 (solid red line), and 7 = 10 (dashed 
blue line). (Top) Positive n. For 7 = 1, models to the right 
of the purple dashed line have large quantum corrections; in 
the 7 = 0.1 case, all models shown pass the quantum sta- 
bility test. For 7 = 10 and n > 3 no models are excluded. 
(Bottom) Negative n. All models shown pass the quantum 
stability test. 



7. For n = 4, the self-coupling 7 = 10 corresponds to 
A = 240, a rather large number for which the chameleon 
effect is very strong. Thus there are no constraints for 
this value in Fig, gg (Top). 

Both plots in Fig. [Tg show that constraints vanish 
at large \n\. We can see why by computing the max- 
imum possible force per unit area between two planar 
slabs of density p m . If the distance separating them is 
small, then the field at the center of the gap will be 
cj) g ss 4>B(Pm)- If each slab is sufficiently thick, then 
the field on the side facing away from the other slab 
will be cj) sv = 0B(/5 m )(f ^ l/ 71 )- Then the magnitude of 
the attractive force between them is F = f3p m \4>B{Pm) — 
<M/Mpi = (/WA/pi)|<Mp m )/H -> /3p m M A /\nM P1 \ at 
large n. Suppose that we also include a force suppression 
factor sech(2m e ff Azfoii). At large |n|, m c s ~ \J 0\n\, so 
the suppression factor decreases quickly. Thus large- |n| 
models will be difficult to exclude. 

As a final note, we have used the lDpp approximation 
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FIG. 13: Field on the surface of a wedge on the source disk, 
for A = 100, f3 = 10, and Az S -t = 0.1 mm. 



to study the effects on these constraints of a degraded lab- 
oratory vacuum. We find that the chameleon fifth force 
at Azs-T = 0.1 mm is extremely insensitive to the vac- 
uum quality; even conducting the experiment at atmo- 
spheric pressure does not noticably reduce the chameleon 
fifth force. At Az s _ T = 10 mm, the largest disk sepa- 
ration probed by Eot-Wash, constraints at atmospheric 
pressure are ~ 10% worse than those in a vacuum for 
n = -1. 



D. Forecasts 

The geometry of the next-generation Eot-Wash source 
and test disks is shown in Figure [6] (Right). In order to 
visualize such a disk, one can imagine a pie cut into 240 
equal wedges, with every other wedge removed, and a 
circular region excised from the center. We approximate 
each disk as having an inner radius of 13 mm and an outer 
radius of 23 mm. Each of the N we d ge s = 120 wedges has 
a thickness Az w = 50 /jm and a height Ay w = 10 mm. 



As in 47] , we approximate each wedge as a rectangular 
sheet of width Ax w = 27rr aV g/(2iV wcc jgcs) = 0.47 mm, 
where r avg = 18 mm is the average of the inner and 
outer radii. Each wedge has a density p m = 20 g/cm 3 , 
and they are mounted on a glass disk of density 2 g/cm 3 . 

Figure [13] shows this lDpp approximation along with 
the more accurate three-dimensional numerical simula- 
tion. Agreement between the two is not as close as it 
was in Figs. [7JIH1 This is because edge effects are larger 
when the features in the disks are long, narrow grooves 
rather than circular holes. However, our approximation 
reproduces the qualitative features of the field, and par- 
ticularly the difference in the surface field between re- 
gions which do and do not overlap a wedge on the op- 
posite disk. This field difference determines the energy 
difference, hence the predicted torque. 

Let the gap between source and test disks be Azs-T- 
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FIG. 14: Torque as a function of rotation angle for the next- 
generation Eot-Wash apparatus, assuming A = f3 = 1 and 
Azs-t = 0.1 mm. The lDpp approximation overestimates 
the numerical computation of [47[ by ~ 50%. 



Wc can immediately apply the lDpp approximation to 
determine the energy associated with the overlap between 
wedges on opposite disks. 



E = J dx dy dz 

_ Ay w f z 
2m, 



V0 



ff Jo 



dx 



2m 



cir \ 



d(f> s 


2" 


dx 





(40) 



where f z = 1 — exp(— m e ffAz w ) corrects for the finite 
wedge thickness; as in (|37]). m c s and 4>b are evaluated 
at the bulk density p m . In the case of perfect overlap, 
4> s (x) is a constant equal to the surface field in a gap of 
size Azg_x- In the case of no overlap, 4> s (x) is the surface 
field in a gap of size yj (Azs_t) 2 + (Ax w /2 — x) 2 . 

After integrating to find the energy difference AE be- 
tween the overlapping and non-overlapping configura- 
tions, we may proceed as before to find the torque, 



r = N wcd gcs fsu P AE sin(iV wedges 6 



(41) 



This lDpp approximation is compared to the three- 
dimensional numerical calculation in Figure 1141 lDpp 
overestimates the correct torque by « 50%. This is likely 
due to Approximation [2] in Sec. IIV A[ which artificially 
flattens out the field in the region x > in Fig. [T3l which 
overlaps the test wedge. This leads to an overestimate of 
the energy difference and torque. 

As with the current experiment, we assume that the 
next-generation experiment places an upper bound on 
the torque of 0.003 fN-m at a separation distance Azg_T = 
0.1mm. We also assume an identical shielding foil. The 
resulting forecasts are shown in Figure [T5] (Top) for the 
</> 4 potential. In particular, we note that for a range of 
matter couplings 0.1 < f3 < 1000, Eot-Wash will be able 
to exclude all <^ 4 chameleon models satisfying the quan- 
tum stability condition of Sec. Ill CI This is an improve- 



large quantum corrections 




Eot-Wash 
next-generation 
(1 Dpp approx.) 



10 100 
matter coupling p 



1000 10000 



10 4 



10 J 



10 



10' 



10 u 



10" 



10"' 



10"' 




Eot-Wash 
next-generation 
(1 Dpp approx.) 



corrections 



0.01 



0.1 



1 10 
matter coupling p 



100 



1000 



FIG. 15: Forecast constraints from the next-generation Eot- 
Wash apparatus. The light green shaded region is ex- 
cluded; models in the black region are linear inside the disks. 
(Top) n = 4. Models inside the long-dashed blue curve are 
excluded by the current Eot-Wash experiment, while models 
above the short-dashed purple line have large quantum cor- 
rections. (Bottom) n = — 1. Models below the short-dashed 
purple curve have large quantum corrections. 



ment of two orders of magnitude relative to the current 
experiment. 

Figure [T5] (Bottom) forecasts constraints on the n = 
— 1 chameleon. Constraints themselves are not substan- 
tially stronger than those of the current experiment. 
However, differences in the geometry and the density 
mean that quantum corrections are larger. All quantum- 
stable, nonlinear n = — 1 chameleons with 0.07 < (3 < 5 
will be excluded by this experiment. 

Constraints on models with large \n\ at 7 = 1, shown 
in Figure [TBI will improve dramatically relative to those 
of the current experiment. In the case of n > 2, 
the next-generation Eot-Wash will exclude all quantum- 
stable chameleon models with 0.1 < /3 < 1000. For 
n < —1, it will exclude all such chameleon models with 
0.1 < P < 20. This is consistent with the claim of [5l| 
that an order-unity improvement relative to the current- 
generation experiment would allow Eot-Wash to exclude 
all quantum-stable chameleon models with matter cou- 
plings near unity. 
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curve have large quantum corrections. (Bottom) Negative n. 
Models below the curve have large quantum corrections. 



Since the lDpp approximation underestimates the 
torque in the current experiment but overestimates it in 
the next-generation experiment, it is possible that the 
forecasts presented here are an overestimate. As dis- 
cussed in the previous section, we could introduce a cor- 
rection factor to be calibrated using the numerical so- 
lutions of [47|. It is also possible that the experimental 
sensitivity is somewhat worse than that of the current ex- 
periment, or that the next-generation Eot-Wash probes 
a somewhat different distance scale Azs_t- Figure [17] 
shows the constraints which would result for 4> 4 theory if 
the distance or sensitivity differ from our assumed values. 



V. CONCLUSION 

Modern torsion pendulum experiments, designed to 
test Newtonian gravity on submillimeter distance scales, 
are capable of uncovering new physics at the dark en- 
ergy scale of M A = 2.4 x 1CT 3 cV ~ (1/0.1 mm). We 
have developed an approximation allowing us to estimate 
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FIG. 17: Dependence of lDpp forecasts on the torque bound 
(previously assumed to be r max = 0.003 fN-m) and the dis- 
tance probed (previously assumed to be Azs-t = 100 /j,m). 



the chameleon-mediated fifth force which would result in 
a torsion pendulum experiment such as Eot-Wash as a 
function of the chameleon model parameters. This is ac- 
complished by mapping the geometry of the source and 
test masses locally onto a one-dimensional plane-parallel 
problem, which can be solved exactly in a series ex- 
pansion. This approximation accurately reproduces the 
chameleon field on the surface of each mass and allows 
us to compute the expected torque signal as a function of 
rotation angle, correct to a factor of ~ 2. Furthermore, it 
agrees well with published constraints on 4 chameleon 
fifth forces using the current-generation Eot-Wash exper- 
iment. 

We have used this approximation to extend Eot- 
Wash constraints to a much wider range of chameleon 
models, as shown in Figures [IT] and [12] Of particu- 
lar interest is the quantum stability condition described 
and Sec. Ill C[ which argues that current 
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torsion pendulum experiments are on the verge of ex- 
cluding all chameleon models with small loop correc- 
tions and gravitation-strength matter couplings j3 ~ 1. 
The current experiment can exclude such quantum-stable 
chameleons for certain specific models, but constraints on 
them remain weak for inverse power law self-interactions. 

Additionally, we have forecast constraints from the 
next-generation Eot-Wash experiment. This experiment 
is expected to be powerful enough to exclude a large 
range of models satisfying the quantum stability condi- 
tion. We show in Fig. [TS] (Top) that the next-generation 
Eot-Wash will exclude all quantum-stable 4 chameleons 
with matter couplings in the range 0.1 < (3 < 1000, 
an improvement by two orders of magnitude relative to 
the current experiment. For unit self-interactions, the 
next-generation experiment will exclude all quantum- 
stable n > 2 models with 0.1 < (3 < 1000 and all such 
n < — 1 models with 0.1 < j3 < 20, as illustrated in 
Fig. [T5] With the potential to detect or exclude a vast 
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range of quantum-stable power-law chameleon models 
with gravitation-strength couplings, the next-generation 
Edt-Wash experiment will be a powerful probe of dark 
energy candidates at the laboratory scale. 
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